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ABSTRACT 

Using archival data from the Hubble Space Telescope, we study the quantitative mor- 
phological evolution of spectroscopically confirmed bright galaxies in the core regions 
of nine clusters ranging in redshift from z = 0.31 to z = 0.84. We use morphological 
parameters derived from two dimensional bulge-disk decomposition to study the evo- 
lution. We find an increase in the mean bulge-to-total luminosity ratio B/T as the 
Universe evolves. We also find a corresponding increase in the fraction of early type 
galaxies and in the mean Sersic index. We discuss these results and their implications 
to physical mechanisms for evolution of galaxy morphology. 
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1 INTRODUCTION 

Clusters serve as laboratories for investigating the depen- 
dence of galaxy morpho logy on the density of the envi- 
ronment. iDressler] l|l980l ) found that the fraction of early 
type galaxies increases with local density of galaxies. He 
also found that around 80% of galaxies in nearby clusters 
are of the early type. In recent years, deep, high resolution 
imaging by the Hubble Space Telescope (HST) has helped 
to extend the study of galaxy morphology as a function of 
the environment to z ~ 1. The observed fraction of differ- 
ent morphological types can be explained in terms of both 
'nature' and 'nurture' scenarios. The former argues that a 
galaxy's morphological typ e is determined by initial con- 
dition s at formation (e.g. lEggen. Lvnden-Bell fc Sandagd 
1962), while the latter depends on the influence of the en- 
vironment and of secular e volution for de termining the final 
morphological type (e.g. iToomre! 1 19771 ). Numerical simu- 
lations play an important role in investigating the relative 
importance of the two scenarios under different physical con- 
ditions. 

Following Dressler's pioneering work, many observa- 
tional studies have been done to measure and understand 
the morphology density relation (MDR), and the depen- 



* E-mail: vvinuv@iucaa.ernet.in 
f E-mail: yogesh@ncra.tifr. res. in 
| E-mail: akk@iucaa.ernet.in 
5 Deceased 



dence of morphologi cal type on distance of the galaxy from 



Goto et al.l 



Holden et al.l 



the cluster centr e iTDressler et al.l Il997 
jSmith et all 120051 ; iPostman et al.ll2005l : I 
iTreu et al.1 l2003t ICapak et all 120071) . ISmith et all (|2005D 

found that the early type fraction is constant in low density 
environments over the last 10 Gyr, but there is significant 
evoluti on in this fraction i n higher density regions. Accord- 
ing to ISmith et al.l l|2005t ). this suggests that most of the 
ellipticals in clusters formed at high redshift, and the in- 
crease in the fraction of early type galaxies is because of 
the physical processes in dense regions which transform disk 
galaxies with ongoing star formation to early types. Dissipa- 
tionless merging of cluster galaxies may also be responsible 
for this increase. 

There is increasing evidence to show that massive el- 
lipticals formed by dissipationless (dry) merger of two or 
more systems. Such ellipicals must have f ormed at later 
times t han their low l umino sity counterparts (Ide Lucia et al.l 
120061 ). Ivan Dokkuml l|2005l ) analysed tidal debris of ellipti- 
cal galaxies and concluded that ~ 70% of the bulge domi- 
nated galaxies have experienced a merger. The analysis of 
nearby bulge dominated galaxies has shown that the gas 
to stellar mass ratio is very small and these mergers are 
mostly 'dry'. Using a semi-analytic mod el for galaxy for- 
mation, iNaab. Khochfar fc Burkertl (|2006l ) found that both 
the photometric and kinematic properties of massive ellipti- 
cal galaxies are in agreement with the scenario where mas- 
sive elliptical galaxies are produced by mergers of lower 
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mass ellipticals. They suggested that the merger of two 
spiral galaxies alone cannot reproduce the observed prop- 
erties, and that the large remnant mass (> 6 x 10 Mq) 
implies that they must have undergone elliptical - ellipti- 
cal mergers. They found that this process is independent of 
the environment and redshift, which means that dry merg- 
ers can occur at low redshift as well. By analysing spirals 
from cluster an d group environmen ts at intermediate red- 
shift (z ~ 0.5), iMoran et all (120071 ) found that the Tully- 
Fisher relation shows larger scatter for cluster spirals than 
for those in the field. They also found that the central sur- 
face mass density of spirals in clusters is small beyond the 
cluster virial radius, and argued that these observations pro- 
vide evidence for merger/harrasment. 

Mergers are not common in clusters, as the velocity 
dispersion of virialized clusters is large. So if ellipticals in 
clusters have formed by mergers, that most likely happened 
durin g the early stage of cluster collapse (|Roos fc Aarseth] 
1982). There is some observati onal evidence to support this 
idea; Ivan Dokkum etail (|l999h showed that there is a large 
fraction of ongoing mergers in a cluster at z — 0.83. The 
observation of the unvirialized cluster RX J0848 + 4453 at 
z = 1.27 revealed many ongoing d issipationless mergers of 
galaxies l|van Dokkum et aL [2001) . These observations go 
against the view of monolithic collapse where all the ellipti- 
cals formed at the same time, at very high redshift. 

In this letter we report on the evolution of galaxies in 
the core region of clusters, measured using bulge-disk de- 
composition. We show that in the case of brightest clus- 
ter galaxies, the fraction of galaxies with B/T > 0.4 and 
n > 2.5 evolved significantly over the redshift range 0.31 
to 0.83. Throughout the paper we use the standard concor- 
dance cosmology with S1a = 0.73, Q. m = 0.27 and Ho — 71 



km : 



Mpc 



2 CLUSTER SAMPLE AND DECOMPOSITION 
TECHNIQUE 

We have exclusively used archival Hubble Space Telescope 
(HST) data in this work. All of the observations used by 
us were obtained with either the ACS or WFPC2 cameras 
on-board the HST. We constructed our sample of clusters 
by an extensive literature survey of HST observations of 
moderate redshift clusters. The clusters we study span the 
redshift range 0.31 to 0.837. The clusters were selected such 
that they each had at least 15 spectros copically confirmed 



cluster members listed in the literature (Couch et al 



Smail et al.lll997l; iDressler et al.|[l999l; IWilliams et al 



1998; 



1999; 



Halli dav et al.ll2004l ; lDemarco et al.ll2005l ). We restricted our 



study only to those clusters whose second brightest member 
galaxy has an absolute magnitude between -25.0 and -27 in 
rest-frame B-band. Here the magnitudes are corrected for 
the cosmological surface brightness dimming. This criterion 
allows us to restrict ourselves to clusters of roughly com- 
parable luminosity. Further, we only included clusters with 
imaging data around the rest-frame B filter i.e. in filters 
where the central wavelength of the filter corresponded to 
a rest frame wavelength in the range 350-550 nm. This has 
the advantage that the k-correction for transformation from 
Rot I filters in the observer frame, to the B filter in the rest 
frame, has only a weak dependence on galaxy spectral type 



in the redshift range 0.3 < z < 0.8 l|Bohm fc Zieglerll2007f ). 
Incompleteness of the spectroscopy is a potential problem; 
since our data are drawn from a heterogeneous set of ob- 
servations, it is not possible to correct for incompleteness in 
a completely consistent way. However, for the five clusters 
in our sample whi ch are part of the ES Q Distant Clusters 
Survey (EDisCS), lHallidav et all l|2004h have shown that 
incompleteness does not introduce a significant bias. The 
level of incompleteness is also relatively small because they 
(like us) restrict themselves to the brighter cluster members. 
Nine clusters satisfied our selection conditions at the time 
of commencing this study; the basic data for these clusters 
are given in Table [TJ 

The images we downloaded from the HST data archive 
were processed in a standard way using the On-the-Fly Re- 
processing (OTFR) pipeline at STScI. The processed im- 
ages were dark and bias subtracted and flat-fielded by the 
OTFR p ipeline. We combined t hem using the Multidrizzle 
package (|Koekemoer et al"1l2002l ) to flag and remove cosmic 
rays and to correct for geometric distortion and produced 
a single coadded image. For some clusters, multiple disjoint 
pointings have been used; in such cases we obtained multiple 
multidrizzled images. The correlated pixel noise introduced 
by the drizzle process, was cor rected for, using the prescrip- 
tion of Cas ertano et alj (|2000h . 

We computed the center of each cluster as the centroid 
of the brightest cluster galaxy. We then selected all galaxies 
with a spectroscopic redshift confirming their cluster mem- 
bership and located within a 1 Mpc projected distance from 
the cluster center. Our study is therefore restricted to galax- 
ies lying (mostly) within the core region of the clusters. Note 
that for three clusters - AC 114, CL 0303+17, and 3C 295 
- the HST imaging is not complete for the 1 Mpc projected 
distance. 

Galaxies are known to undergo luminosity evolution, to 
account for which we adopted the following simple scheme: 
M v (z) = M v (z = 0) - 0.8z, wh ere z is the redshift of the 
object and M v (z = 0) = -19.5 l|Postman et alj|2005h . The 
magnitude cutoff in the observed HST filter was then calcu- 
lated using 



m Um = M v {z) + DM - (My - Must) + k 



'HST 



(1) 



where DM is the distance modulus, Mv — Mhst is the rest- 
frame color between V and the observed HST filter and 
kn st is the k-correction in the observed HST filter. The 
k-correction for each cluster was comp uted using the ellip- 
tical SED provided bv lPoggiantil (|1997m . All magnitudes are 
in the Vega system. 

With all the above constraints, we obtained a sam- 
ple of 379 galaxie s in nine clusters. We used the GALFIT 
[|Peng et alj|2002h program for 2D bulge disk decomposition 
of the galaxy images. To fit the galaxies, we have developed 
an automated pipeline (Vinu et al. 2010, in preparation) 
which completely automates the fitting procedure and orga- 
nizes the results into a database, supplemented by a number 
of diagnostic plots to detect anomalies. The basic steps in 
performing the 2D decomposition are: (1) make a cut-out 
image for each galaxy, including neighbouring galaxies. (2) 
estimate starting values of fitting parameters using the pa- 
rameters measured by Sextractor (3) make masks to exclude 
stars, faint galaxies and other artifacts and (4) Run GAL- 
FIT. Bright neighboring galaxies were fitted simultaneously 
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with the target galaxy. We modeled each galaxy as a linear 
sum of a Sersic (for the bulge) and an exponential (for the 
disk) component. The centers, position angles, ellipticities 
and central surface brightnesses of bulge and disk were fit- 
ted simultaneously. The Sersic index was left unconstrained. 
The estimation of galaxy parameters, particularly the Sersic 
index, may be systematically affected if the sky or PSF are 
incorrectly estimated. To test for the effect of using an in- 
correct PSF, we ran the decomposition for each galaxy using 
two stellar PSFs, one constructed using the nearest star and 
the other using the second nearest star and compared the 
fitted values of B/T and n for every galaxy. To test for in- 
correct sky, we ran the decomposition in two modes: one in 
which the sky was left as a free parameter and one in which 
it was fixed to the local sky value as determined by Sextrac- 
tor. In both tests, ~ 80% of galaxies showed random changes 
in the extracted parameters at < 10% level, with no obvious 
systematics. Our s imple tests are consis tent with the exten- 
sive simulations of lHaussler et all (|2007h which showed that 
GALFIT estimates parameters accurately for HST data, 
even at relatively shallow depths. As an additional precau- 
tion, we examined the fit results, by eye, for every galaxy. 
The fit diagnostic plots were used to evaluate the quality of 
the fit. We used the reduced \ 2 given by GALFIT to identify 
galaxies with large residuals. A large residual in the central 
part of the galaxy may be caused by improper estimation of 
the point spread function (PSF). Also, simultaneous fitting 
with neighbor galaxies needs special care as the number of 
free parameters is significantly larger. For galaxies that show 
a large residual, we took these caveats into consideration and 
refit the galaxies until the residual became small. We then 
assigned a quality factor for each galaxy based on the mag- 
nitude of the residual and the deviation of fitted position 
angle from the galaxy position angle estimated by eye. If a 
fit was below a threshold quality, we excluded that galaxy 
from further analysis. We found that most of the galaxies 
that failed the quality check, either have peculiar morphol- 
ogy or strong spiral arms. We successfully fit 337 out of 379 
galaxies in our sample. Unless otherwise stated, all further 
discussion in this letter only applies to these 337 galaxies. 
The cluster-wise breakup of galaxies with a good fit is given 
in Tabled] 



3 RESULTS 

3.1 Evolution of mean bulge-to-total luminosity 
ratio {(B/T)) and Sersic index n 

In Table [2] we list mean and median values of a few pa- 
rameters of interest for all the nine clusters. Note that the 
bulge-to-total luminosity ratio is computed using the pa- 
rameters of the best fit model. We find the (B /T) of cluster 
galaxies in the central 1 Mpc of the clusters changes from 
0.59±o;o2 at redshift z = 0.31 to 0.48 ± 0.03 at z = 0.837 
(Figure[T]) . Err ors were measured using th e bootstrap resam- 
pling method (lEfron fc Tibishiranil [l993l ) . in this and sub- 
sequent figures. The increase in (B/T) ratio has been found 
qualitatively (i.e. u sing visual morphological classification) 



by previous studie s l|Dressler et aljfl99Tl Fasano et al.ll200C 
ISmith et al.ll2005T ); the present work obtains the result quan- 
titatively using bulge-disk decomposition. It must be noted 



Table 2. Mean and median parameters obtained through bulge- 
disk decomposition 



Cluster 



(B/T) 



(n) 



n 



AC 114 
CL 0303 
3C 295 
CL 1232-1250 
CL 1054-1146 
CL 1040-1155 
CL 1054-1245 
CL 1216-1201 
RX J0152.7-1357 



0.59 
0.52 
0.52 
0.48 
0.45 
0.42 
0.47 
0.48 
0.48 



+ 0.03 

0.02 
+0.05 
0.05 
+0.04 
0.05 
+0.04 
0.04 
+0.05 
0.05 
+ 0.08 
0.07 
+0.05 
0.05 
+0.03 
0.04 
+0.03 
0.03 



4.13 1 



,+0.43 
-0.39 

4 24 +i.02 



■AAV 



3.47 
4.15 
3.60 



0.79 
+ 0.5!) 

0.56 
+0.56 
0.52 
+0.76 



7+0.27 
-0.30 

3.31+ - 20 



-0.65 

3.50+ 06 " 



2 
3.25 
2.74 



0.68 

+0.38 
-0.41 
+0.54 

0.44 
+0.28 

0.30 



3.17 
3.30 
2.41 
3.25 
2.55 
2.57 
2.49 



0.83 
+0.72 
0.78 
+0.41 
0.52 
+1.02 
0.43 
+0.71 
0.56 
+0.56 
0.55 
+0.40 
0.51 
+0.43 
0.28 



0.55 

0.46 
0.43 
0.47 
0.26 
0.28 
0.37 
0.32 
0.40 



+0.03 
0.04 
+0.05 
0.07 
+0.03 
0.04 
+0.04 
0.05 
+0.05 
0.06 
+ 0.08 
0.08 
+0.05 
0.06 
+0.04 
0.04 
+0.02 
0.02 



(B/T) : mean value of B/T bulge-to-total luminosity ratio, (n) 
mean value of Sersic index, n : median value of Sersic index, /;, 
fraction of bulge-like galaxies (see text for definition) 



0.37 0.43 



0.74 0.83 




4.0 4.5 5.0 5.5 6.0 6.5 7.0 7.5 
Lookback Time (Gyr) 

Figure 1. Evolution of mean value of the bulge-to-total luminos- 
ity ratio B/T 



that, if our sample is somewhat biased towards luminous red 
galaxies at high-z, and if such galaxies are early type even at 
those redshifts, then our estimated (B/T) value represents 
an upper limit at z ~ 0.8. In addition, incompleteness will 
also tend to make our estimate of (B/T) too high because 
at fainter magnitu des we would prefe rentially miss low B/T 
late- type galaxies (|Desai et alj|2007l ). This leads to our es- 
timated (B/T) to be the upper limit for these clusters. So 
the selection bias, if any, will lead to an apparent weaker 
evolution. So our estimation of evolution of (B/T) may be 
an underestimate. 

We also find that the mean value of Sersic index de- 
creases with lookback time over this redshift range. The 
value changes from 4.13± jg to 2.74j;° ^ from z = 0.31 
to z = 0.84. Figure [2] shows the evolution of Sersic index 
against lookback time. 



3.2 Evolution of bulge-like galaxy fraction 

Since the (B /T) decreases significantly with lookback time, 
one expects to see a simultaneous decrease in the fraction of 
early-type galaxies with a bulge-like morphology. We clas- 
sify a galaxy as bulge- like if its B/T ^ 0.4 and n > 2.5. The 
second condition is required to exclude galaxies with a very 
strong disk which can, on occasion, be incorrectly modeled 
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Table 1. Summary data for sample clusters 



Cluster 


RA 


Dec 


z 


Camera 


Filter 




N T ot 




Reference 


AC 114 


22 58 48.4 


-34 48 60 


0.31 


WFPC2 


F702W 


20.77 


72 


68 


Couch et al. (1998) 


CL 0303+17 


03 06 15.9 


+17 19 17 


0.42 


WFPC2 


F702W 


21.64 


28 


26 


Smail et al. (19971 


3C 295 


14 11 19.5 


+52 12 21 


0.46 


WFPC2 


F702W 


21.93 


37 


32 


Smail et al. (19971 


CL 1232-1250 


12 32 30.3 


-12 50 36 


0.54 


ACS 


F814W 


21.66 


16 


41 


White et al. (20051 


CL 1054-1146 


10 54 24.4 


-11 46 19 


0.697 


ACS 


F814W 


22.44 


30 


25 


White et al. (20051 


CL 1040-1155 


10 40 40.3 


-11 56 04 


0.704 


ACS 


F814W 


22.47 


25 


19 


White et al. (20051 


CL 1054-1245 


10 54 43.5 


-12 45 51 


0.75 


ACS 


F814W 


22.68 


29 


28 


White et al. (20051 


CL 1216-1201 


12 16 45.3 


-12 01 17 


0.794 


ACS 


F814W 


22.89 


50 


43 


White et al. (20051 


RX J0152. 7-1357 


01 52 27.4 


-13 55 01 


0.837 


ACS 


F775W 


23.50 


62 


55 


Blakeslee et al. (2006") 


miim '■ Faint maj 


;nitude cutoff 


in the observed HST filter, N Tot 


: Total 


number 


of galaxies, Npu 


: Number of 



galaxies with good fit 



0.31 0.37 0.43 0.50 0.57 0.65 0.74 0.83 0.94 




3.5 4.0 4.5 5.0 5.5 6.0 6.5 7.0 7.5 
Lookback Time (Gyr) 

Figure 2. Evolution of mean value of Sersic index 



as a bulge with n ~ 1 (an exponential), with a correspond- 
ingly high B /T luminosity ratio. With this definition, we are 
able to compute a bulge-like galaxy fraction for each clus- 
ter. In figure Owe plot this fraction against lookback time. 
Note that the fraction is normalized by the total number 
of galaxies in the cluster, not the number of galaxies with 
a successful fit. The galaxies that are poorly fit are dom- 
inated by galaxies of irregular morphology. We see a near 
monotonic decrease with lookback time in the bulge domi- 
nated fraction of galaxies. We find 40.0 t% % °f galaxies at 



redshift z = 0.837 are bulge-like. This increases to 55 
within ~ 3.5 Gyr. 



+3 



% 



3.3 Discussion 

In the last decade, several studies have reported evolution 
of the morphological content of galaxy clusters by a visual 
study of galaxy morphology (|Smail et all 19971 ; I Couch et al.l 
119981 ; iFasano et alj I2OO0I ; iDesai et alj|2007h . Other studies 



focused on the evolution of the morphology-density rela- 
tion (MDR), where morphology chang es were studied as a 
function of redshift and galaxy density dDressler et al.ll 19971 ; 
iTreu et al.ll2003l ; IPostman etalll200a ISmith et alj|2005h . 

In this work, we have taken the first approach; however, 
we have used quantitative measures of galaxy morphology 
rather than a qualitative morphological classification by eye. 
We decompose the galaxy light into bulge and disk compo- 
nents and study the evolution of morphology of galaxies in 
clusters. We find that the bulge component of the bright 



galaxies in clusters is, on average, becoming stronger as the 
Universe evolves. The fraction of bulge-like galaxies, defined 
as having B/T > 0.4 and n > 2.5 also increases from 40% to 
55%. There is also a significant increase in the mean value 
of the Sersic index as the Universe ages. 

Our results on the evolution of morphological frac- 
tion are consistent wi t h previous work on the sub- 
ject dDressler et all 1 19971; IPostman et al.ll2005l ; ISmith et all 
1200511. However, IDesai et alj l|2007l lT using ACS observations 
of galaxy clusters with 0.5 < 2 < 0.8 found no evolution of 
morphological fraction. At first glance, this seems to con- 
tradict our results. However, it must be noted that in our 
study, we are only including bright, spectroscopically con- 
firmed cluster galaxies. If fainter galax i es are included, as 
was done in the sample of IDesai et al.l (|2007t ) . the ~ 10% 
change we see in the morphological fraction may easily be 
washed out. This ex planation also agrees with the results of 
iHolden et all l|2007f ) who showed that there is evolution in 
the early type fraction with redshift if a luminosity limited 
sample is used. 

IPostman et alj (|2005T l; IDesai et al.l l|2007l l; 
iPoggianti et al.l l|2009l l found that the E+S0 fraction 
correlates with the velocity dispersion of the cluster. The 
fraction is large for clusters with high dispersion. This effect 
would be consistent with the evolution we see, provided 
velocity dispersion systematically decreases with redshift. 
We see no systematic dependence of velocity dispersion 
with redshift except that the two highest redshift clusters 
in our sample have high velocity dispersions (> 1000 km/s). 
This high value, which implies higher (B/T) is consistent 
with the mild increase of (B/T) and the fraction of bulge 
dominated galaxies for the two high 2 clusters. But the 
dominating effect seems to be morphological evolution. 

The (B/T) of a cluster may increase either due to an 
increase in the strength of the bulge component of galax- 
ies or by the fading of disk component. The increase in 
average bulge strength may be caused by me rging which 
tends to produce a elliptical like morphology ijvan Albadal 
1 19821 ; lBarnes!ll992l ; lHernquist|[l992l ; iBournaud et al.ll2005T l. 
On the other hand, the fading of the disk could be a 
by-product of the morphological transformation of galax- 



In clust ers, a variety of mechanisms such as 



I Moore et al 



harassment ( Moore et al. 1 19961 ; iMoore. Lake fc Kat 



alax 



199 



19991). minor mergers and ram pressure strip- 



ping (|Gunn fc Gottl 1 19721 ; lAbadi. Moore fc Bowerl I1999T ) 



may contribute to the disk fading. 
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0.31 0.37 0.43 0.50 0.57 0.65 0.74 0.83 0.94 




3.5 4.0 4.5 5.0 5.5 6.0 6.5 7.0 7.5 
Lookback Time (Gyr) 



Figure 3. Evolution of fraction of bulge-dominated galaxies 



The increase in the mean Sersic index (n) with cos- 
mic time seems to indicat e that mergers play a role 
l|Scannapieco fe Tisseral l 200St). Merger event s are common 
at intermediate redshift dDressler et al.|[l994l ). The increase 
in the fraction of galaxies with high B/T and Sersic index 
is possible if galaxies gradually evolve into a phase where 
the spheroidal component increasingly dominates. The end 
point of such evolution is an elliptical galaxy. 

It has recently been suggested that the bright SO popu- 
lation has lik ely formed through monol ithic collapse or ma- 
jor mergers l|Barwav et al.l [20071 . [2009). Numerical simula- 
tions have also shown tha t dissipative merger of two un - 
equal mass disk galaxies i|Bekkif Il99sl ; iBekki et all 120051 ) 
can produce lenticulars. Recent observations of low redshift 
clusters suggest that a majority of the infall population is 
merging or interacting ([Moss 2006). Coupled with the con- 
straints from ICarlbera (|l986t ) , it is becoming increasingly 
clear that formation of spheroidal (mostly lenticular) galax- 
ies through mergers is the dominant mechanism behind sys- 
tematic changes in galaxy morphology. 

It must be noted that the trends we see are weak and 
statistical in nature; they are only visible when averaged 
over a large number of galaxies in a large number of clusters 
over a wide range of redshift. The detailed physics operating 
in each cluster, doubtless modifies the morphological evolu- 
tion of galaxies in that cluster. Nevertheless, the fact that 
we see trends indicates that they are real and strong enough 
not to be drowned by the different physical conditions and 
processes operating in individual clusters. Using the large 
database of bulge/disk decomposition results we have ob- 
tained, we are attempting to disentangle cluster specific ef- 
fects from cosmological ones. 
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